Inhomogeneous superconducting phases in the frustrated Kondo-Heisenberg chain 
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We use bosonization and renormalization group methods to determine the ground state phase 
diagram of a one-dimensional frustrated Kondo-Heisenberg system consisting of a one-dimensional 
spin-1/2 Luttinger liquid coupled by a Kondo exchange interaction Jk to a frustrated quantum an- 
tiferromagnetic Heisenberg chain, with a nearest-neighbor exchange coupling J\ and a next-nearest- 
neighbor (frustrating) exchange interaction J2 . We analyze the interplay of quantum frustration in 
the antiferromagnetic chain with the Kondo exchange coupling Jk with the Luttinger liquid. We 
discuss the structure of the phase diagram of this system as a function of the ratios Jk / J\, J2/ Ji 
and of the parameters of the Luttinger liquid. In particular we discuss in detail the regimes in which 
a pair-density-wave state may be realized and its relation with the spin correlations in the frustrated 
antiferromagnetic chain. 
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I. INTRODUCTION 

Kondo-Heisenberg chains arc simple model systems in 
which a (one-dimensional) Luttinger liquid is coupled to 
a one-dimensional antiferromagnetic system by a Kondo 
exchange interaction. Systems of this type display a 
variety of non-trivial ground states with a dazzling ar- 
ray of unconventional behaviorsji^~ In particular, they 
have phases that exhibit a spin gap and have supercon- 
ducting correlations with the peculiar aspect that their 
only viable order parameters are made of observables 
from the spin chain and of the Luttinger liquid which 
separately have only short range correlations. Kondo- 
Heisenberg systems constitute an ideal testing ground in 
which the physics of strongly correlated electronic sys- 
tems can be studied using controlled analytic approxi- 
mations and powerful numerical methods. 

One of the phases that is well established to exist in 
the Kondo-Heisenberg chain is a state with a spin gap 
and long-range (power-law) superconducting correlations 
with finite momentum (which is commensurate with the 
lattice of Kondo-Heisenberg spins). This state was iden- 
tified in Ref. [7| as a pair-density- wave superconducting 
(PDW) state, i.e. a Larkin-Ovchinnikov state without 
an external Zeeman (magnetic) field. In this paper we 
will consider the effects of next-nearest-neighbor antifer- 
romagnetic exchange interactions in the spin chain. In 
the presence of such interactions the spin chain is frus- 
trated and can be regarded as a zig-zag spin ladder. The 
purpose of this paper is to determine the phase diagram 
of this system as a function of the ratio of the Kondo ex- 
change interaction Jk to the nearest-neighbor exchange 
interaction Ji in the spin chain and of the ratio be- 



tween the next-nearest-neighbor exchange interaction of 
the spin chain J2/J1, as well as the interaction coupling 
constants of the Luttinger liquid. 

One motivation for considering the role of frustration 
is to examine the mechanisms that determine the order- 
ing wave- vector Q of a pair-density-wave state. In a one- 
dimensional fermionic system with at least two bands (or 
species of fermions) with attractive interactions it is pos- 
sible to have a Larkin-Ovchinnikov state with any wave 
vector provided that the Fermi points of the two species 
are different. This is possible since in ID there is always 
nesting of the Fermi points. However, if the strongly cor- 
related system has only repulsive interactions this prob- 
lem is more subtle. For example. Berg and coworkers! 
found that in Kondo-Heisenberg chain the ordering wave- 
vector Q of the PDW is the same as the wave-vector of 
the antiferromagnetic order, which in turn is determined 
by the spacing of the Kondo spins. This is true even 
though in such a spin-gap phase the magnetic order is 
short ranged. A similar result was obtained in Ref. [8| 
in the context of doped antiferromagnetic two-leg spin 
ladders, with the only difference that the ordering wave 
vector is determined by spontaneous-symmetry-breaking 
of translation invariance rather than the explicit breaking 
of translation invariance of the Kondo-Heisenberg chain. 
Thus, in both cases the PDW order is commensurate. 

It is apparent that this commensurate state is the re- 
sult of the magnetic origin of the mechanism of the PDW 
order, provided the SU(2) magnetic symmetry is exact. 
This is true even in the case of the frustrated chain which, 
for J2 larger than a critical value, at the classical level 
the frustrated Heisenberg chain has a ground state with 
incommensurate spiral magnetic order. On the other 



hand, at the quantum level, in the case of a chain of 
spin- 1/2 degrees of freedom the ground state is dimerized 
and has short-ranged spin order.— The short range mag- 
netic order becomes incommensurate order for J2 > Ji/2, 
the Majumdar-Ghosh point^, where the dimerization is 
strongest. 

This behavior naturally poses the question of whether 
in a strictly one-dimensional system an incommensurate 
phase is possible at all. A ground stye with incommen- 
surate magnetic order requires the existence of an addi- 
tional gapless collective mode. In a one-dimensional sys- 
tem this is possible only if the effective low energy theory 
has an exactly marginal operator. However, the effective 
field theories of models with an exact SU(2) symmetry 
(and their generalizations) are non-linear sigma models 
with compact target space manifolds, these target man- 
ifolds do not have flat directions and hence their ground 
states can only have massive (gapped) excitations.^- The 
only alternative to this gapped ground state is that the 
generalized non-linear sigma models may have a topolog- 
ical term. In this case the RG fiows drives the system 
to a non-trivial finite conformally invariant fixed point. 
The only available fixed points with SU(2) symmetry are 
Wess-Zumino-Witten (WZW) models with current alge- 
bras SU(2)i; (where k is an integer) (or generalizations 
of these fixed points). However, fixed points of this class 
also do not have any exactly marginal operators. These 
arguments suggest that it is not possible to have a state 
with incommensurate magnetic order with exact SU(2) 
spin invariance in one spatial dimension. The only way to 
circumvent this problem is either to break the magnetic 
symmetry explicitly, from SU(2) down to a U(l) sub- 
group, as in the presence of magnetic anisotropy, which 
allows for magnetic spiral phases fii or to consider a quasi- 
one-dimensional system in which the symmetry-breaking 
is spontaneous in the higher-dimensional system (or by 
an explicit symmetry-breaking by an uniform magnetic 
field)i2i"— Nevertheless, although magnetic frustration in 
an SU(2) invariant system does not lead to an incommen- 
surate PDW phase in the ID case we are discussing here, 
we will see that it does lead to non-trivial and interesting 
phases that will be discussed in some detail in this paper. 

In this paper we investigate the nature of the ground 
state of a frustrated Kondo-Heisenberg model and the 
resulting phase diagram (see Fig. [T]). Our results are 
based on a combination of bosonization methods (both 
abelian and non-abelian) with perturbative renormaliza- 
tion group arguments. Our strategy is to investigate the 
behavior of the system in limiting regimes of the parame- 
ter space. Thus we will first consider the extension of the 
results of Ref. Q account for the effects of weak next- 
nearest neighbor interactions J2 in the spin chain. Next 
we consider the opposite regime in which the nearest- 
neighbor interactions Ji of the spin chain are weak. This 
regime turn out to be quite rich. Finally we consider the 
regime in which Ji is weak compared with the Kondo 
exchange coupling. In this regime we find a finite non- 
trivial fixed point of the Toulouse type, investigated ear- 



lier on for a somewhat different system by Azaria and 
Lecheminant J^ By considering the leading perturbations 
around the Toulouse fixed point we show that it controls 
a stable phase with the character of a fractionalized spin 
liquid. At some finite value of Ji this phase becomes un- 
stable and has a phase transition to a pair-density-wave 
phase. 

This paper if organized as follows. In Section |TT] we in- 
troduce the frustrated Kondo-Heisenberg model and dis- 
cuss how the order parameters of the different phases 
of interest are realized. In Section IIIII we consider the 
weak frustration regime, J2 <^ Ji, which is treated us- 
ing bosonization methods (both abelian and non-abelian) 
and renormalization group calculations. Here we discuss 
the interplay between the PDW phase of the unfrustrated 
system and the dimerized phase. In Section IIVI we dis- 
cuss the opposite regime, J2 3> Ji which can also be 
treated by bosonization and renormalization group meth- 
ods. In this regime the phase diagram turns out to be 
quite rich and the construction of the order parameters 
is non-trivial. Here we find a stable phase with a gap- 
less fractionalized fluid in the spin sector, a fractionalized 
spin liquid phase, and discuss the quantum phase tran- 
sition from this phase to the PDW phase of the weakly 
frustrated regime. We close this paper by discussing the 
resulting phase diagram of this system, shown in Fig. 
[I] in Section |V] as well as several open questions. The 
Toulouse point solution of the strong coupling regime of 
this system is summarized in Appendix [Xj 



II. THE MODEL 

The frustrated Kondo-Heisenberg chains a model of a 
one-dimensional interacting system of spin- 1/2 fermions 
coupled by a Kondo exchange coupling to a one- 
dimensional array of localized spin- 1/2 degrees of free- 
dom whose Hamiltonian is a Heisenberg antiferromagnet 
with both nearest- and next-nearest neighbor exchange 
interactions. When the latter interactions are larger than 
a critical value the resulting magnetic chain is frustrated 
and can equivalently be depicted as a zig-zag ladder. 

The model is described by the Hamiltonian of Eq. (12. ip 
which is a sum of three terms: a) Hp that represents the 
interacting system of fermions, Ea. (|2.2p . b) Hneis for 
the frustrated Heisenberg, Ea. (|2.3|) , and c) Hk that de- 
scribes the Kondo exchange interaction between the spin 
degrees of freedom in these two subsystems, Ea. (l2.4p . 

H = Hp + Hueis + Hk (2.1) 

Hf = -t 5](clc.+i^ + cl+i.c-) + i/int (2.2) 

HhcIs — Jl 2_^ ^« ■ ^«+l + ^2 2_^ ^i ■ ^«+2 (2-3) 

i i.(T 

Hk - Jk J2 S» • «' (2-4) 

i 

In Eq. (|2.2p TJint represent local interactions in the system 



of ID fermions, e.g. a Hubbard interaction U, a nearest- 
neighbor (Coulomb) repulsion V, and a nearest-neighbor 
(Heisenberg) exchange interaction J. Here we will as- 
sume that the ID system of fermions is gapless and hence 
described by a Luttinger liquid with charge and spin ex- 
citations. Thus, the ID fermionic system has a U(l) x 
SU(2) global symmetry (accounting for the charge and 
spin sectors). The Heisenberg spin chain has an SU(2) 
global symmetry. The Kondo coupling between the two 
systems, represented by Hk in Eq. (|2.4p with coupling 
constant Jk, reduces the symmetry to a U(l) x SU(2) 
global symmetry. In principle the lattice spacing of the 
magnetic chain is different than that of the ID system 
of mobile fermions. For simplicity here we will take the 
lattice spacings of the two subsystems to be the same. 
In this case there is no explicit breaking of translation 
invariance. 

In Eq. (|2.2p we introduced a set of fermion creation 
and annihilation operators, cl ^ and Ci^a, at each site i 
for electrons with spin a =t, i- The three components of 
the spin operator of the itinerant electrons at the ith site 
are given by sf 



kT,a.f3 4,a<,|3'^^,P^ where o = 1,2,3. 



Here tr^ a is the (a, /3) element of the ath Pauli matrix. 
Similarly, the operators S^ represent the spins of the frus- 
trated Heisenberg spin chain with coupling constants Ji 
(for the nearest-neighbor exchange coupling) and J2 for 
the next-nearest-neighbor exchange coupling. 

As we noted before the frustrated chain can also be 
regarded as a zig-zag two-leg spin ladder with J2 repre- 
senting the interactions on the rungs. In this picture, the 
frustrated chain is a ID version of an asymmetric trian- 
gular lattice. This picture is appropriate for quasi-one- 
dimensional systems such as Cs2CuCl4.^^ Both couphng 
constants are taken to be antiferromagnetic. Hence the 
spin chain is frustrated. On the other hand, in what fol- 
lows we will take the Kondo coupling to be the same for 
all sites which is appropriate for the frustrated chain. 

Several limits of this model have been considered be- 
fore. In the case in which the Kondo coupling is zero, 
Jk = 0, it is knowE^il that the frustrated chain has a 
quantum phase transition at a critical value J2C of J2. 
For J2 < J2C, frustration has essentially not effect on 
the low energy properties. In this regime, next-nearest- 
neighbor interaction is marginally irrelevant and the sys- 
tem behaves effectively as the spin- 1/2 Heisenberg model 
chain. Thus, in this regime the decoupled spin chain is 
gapless, and hence critical, and exhibits power law cor- 
relations. When J2 exceeds a critical value of the next- 
nearest-neighbor coupling J2C, the interaction becomes 
marginally relevant. The resulting ground state is a spin 
singlet and there is a finite energy gap in the spectrumj^ 
In this regime the translational order is spontaneously 
broken by the appearance of a magnetic dimerization of 
the chain. The order parameter for the dimerized phase 
is 

e = i-iy (S, • S,+i - S,+i • S,+2) (2.5) 

which acquires a nonzero, position-independent value in 



this phase. The value ^ = 0.241 is known from nu- 
merical studies. At ^ = 0.5, known as a Majundar- 
Gosh point, the ground state could be analytically ob- 
tained. It corresponds to the formation of spin singlets 
between each site and one of its neighbors. In this limit 
the dimerization (and the frustration) is largest. Inside 
the dimerized phase there is short-range incommensurate 
spiral-spin order. Since there is a finite spin gap, not only 
electron tunneling between the two chains is suppressed, 
but also the exchange Kondo coupling becomes irrele- 
vant up to a finite value of the order of the spin gap. 
Hence in this regime the gapless ID electron system and 
the frustrated spin chain are effectively decoupled at low 
energies. 

On the other hand, for J2 = and Jk 7^ 0, the system 
is the Kondo-Heisenberg chain. For small Jk the Kondo 
coupling is a marginally relevant perturbation and fiows 
under the RG to a strong coupling fixed point with a 
finite spin gap and short-range commensurate magnetic 
correlationsi^ In this system only composite order param- 
eters have power-law correlations. ^'^'^ It has been realized 
recently^ that in the spin gap phase (usually referred to 
as the "Kondo singlet" regime in the heavy fermion lit- 
erature) the strongest correlations (i.e. with the smallest 
exponent) describe a pair-density-wave state, a super- 
conducting order with finite wave-vector. In the case in 
which the two lattice spacings are the same the ordering 
wave vector is Q — tt. 

Finally, the limit in which Ji = but with J2 and Jk 
finite has also been studied in some detail. This limit 
is quite rich. It has an unstable fixed point at Jk — 0. 
For finite Jk this system fiows to a non-trivial finite in- 
frared stable fixed point which in some ways resembles 
the physics of the multi-channel Kondo problemJ^ii^iiS 
This fixed point is equivalent to two decoupled chirally 
stabilized systems.^*' In the phase governed by this non- 
trivial fixed point the only allowed order parameters (i.e. 
operators with power-law correlations) are also compos- 
ite. 

In the subsequent sections we will examine the phase 
diagram of the full system by expanding about these two 
limits, J2 ^ Ji and J2 3> Ji, with Jk finite (and with 
the parameters of the ID gapless electronic system fixed). 



III. J2 « Ji 

In this section we will look at the effects of the next- 
nearest-neighbor exchange coupling J2 on the physics 
Kondo-Heisenberg chain. In the regime with small Jk 
and small J2 one can look at the effective field theory 
in a naive continuum limit and analyze the role of var- 
ious operators. Bosonization (abelian and non-abelian) 
are very useful tools to understand the physics of this 
regime.i^i^ The low energy excitations of the electronic 
chain can be taken into account by linearizing the free 
fermion band around the Fermi level. In this regime the 
fermionic operators can be written in terms of the con- 



tinuum right and left fields ip^ji{x) and V'ctl(^) (here the 
label 1 denotes the ID system of mobile electrons) as : 



^--e^'"><,fl,i(a;) + e-''=-" 



V'.L,i(x) (3.1) 



where a is the lattice spacing and kp is the Fermi momen- 
tum of the chain. The local spin operator of the electrons 
Sn = \c\^a'^apCni3, Can also bc dccomposcd into a slowly 
varying piece (with Fourier components near zero wave 
vector) and a rapidly varying piece (with Fourier compo- 
nents near 2k p) which represents the spin-density wave 
order parameter. More explicitly 

s„ ^ a [Jl,i(^) + Jfl,i(^) + e-^2fe^"Ni(x) + h.c.|3.2) 
where 



Jl,1 = -ji^lL,l^c,l3^liL,l 



(3.3) 



are the chiral spin currents of the right and left moving 
electrons, and 



1 



,t 



Ni = ■:^'^lR^i(Tcp'ippL,i 



(3.4) 



For general kp the order parameter operator Ni of the 
ID electronic system is a three-component complex vec- 
tor and, for general filling of the ID electronic system, the 
spin-density wave has an ordering wave- vector Q = 2k p. 
In the special case in which the ID electronic chain is 
half-filled, where there is a Mott charge gap, the Kondo- 
Heisenberg chain is a Kondo insulator. In this case the 
ordering wave vector of the ID electronic system is Q = tt 
and Ni is real (self-adjoint) and is a Neel order param- 
eter. We will not discuss this interesting case here and 
we will focus on the case in which the electronic system 
is metallic. 

On the other hand, the degrees of freedom of the 
Heiscnberg spin chain can also be decomposed into slowly 
and rapidly varying components.— In this case the spin 
operators of the Heisenberg chain could also be repre- 
sented by fermionic fields starting from a half filled Hub- 
bard model (with kp = ir) and a gapped (and hence 
frozen) charge degrees of freedom (due to the Mott gap) . 
The decomposition of the local spin operators of the 
Heisenberg chains is: 

Sn-a [Jl,2{x) + Jfl,2(.^) + (-l)"N2(a;)] (3.5) 

where the label 2 now denotes the spin chain. The ex- 
pressions for the chiral spin currents of the Heisenberg 
chain 3l.2 and 3r^2 are the same as the previous ones by 
changing 1 for 2. The (real) Neel order parameter (the 
staggered magnetization) is: 



N2 = i 



V'1k,2'^":/3V'/3L,2 + V'lL,2'^a,/3V'/3fl,2 (3.6) 



The effective Hamiltonian for the low energy regime 
can be determined from the Hamiltonian of Eqs. 



(12. 3p . and (|2.4p . As it is common in ID problems, here 
too in the low energy regime there is a separation be- 
tween charge and spin degrees of freedom and the effec- 
tive Hamiltonian is a sum of two terms, one for the charge 
sector and one for the spin sector, (for a general discus- 
sion on bosonization see, e.g., Ref. [22|.) The charge 
sector is described by a conventional (compactified) bo- 
son (/)c whose Hamiltonian density, He, is parametrized 
by a charge Luttinger parameter Kc and a charge veloc- 
ity Vc, both of which depend in a non- universal way on 
the microscopic parameters of Eas. (|2.2p . (|2.3|) . and (|2.4p 



1 



nc = ^K + 2^c«c(5,0c) 



(3.7) 



where He is momentum canonically conjugate to the field 
(jjc- The (normal-ordered) charge density jo and charge 
current are related to the field (pc by the usual bosoniza- 
tion formula. 

The Hamiltonian density for the spin sector, "Hs, can be 
written in terms of the right and left moving spin currents 
of the magnetic and electronic chains.^ We obtain 

T~Ls = %(] + T-Lint 

"Ho = — - — (: Jfl,i • J_R,i+: Jl,i • Jl,i : ) 

H -z — (: J7?„2 • J_R,2: +: Jl,2 • Jl,2: ) 

T-Lint = .91 Jfl,! • Jl,1 + 52Ji<;,2 ' Jl,2 
+ .93(Ji?,l • Jl,2 + Jl,1 • Ji?,2) 

-f .g4(Jfl,i-Jfl,2+JL,i-JL,2) (3.8) 

where 1 and 2 label the electronic system and the spin 
chain, respectively, and vi ~ 2at and V2 ^ aJi. All the 
operators that we have included are marginal, and they 
are marginally relevant for g > Q and marginally irrele- 
vant otherwise. Notice that in the effective interaction 
of Eq. ()3.8p we have not included a possible coupling 
between the spin-density wave order parameters of the 
ID electronic system and of the Heisenberg spin chain. 
For general filling of the ID electronic system this inter- 
action is not allowed (or, rather it is strongly irrelevant) 
since the two systems have different ordering wave vec- 
tors. However it is a relevant perturbation in the case of 
the Kondo insulator where it plays a key role. 

The coupling constant g2 parametrizes the strength of 
the backscattering term of the magnetic chain. Accord- 
ing to the discussion of Section [H] this coupling is irrel- 
evant in the absence of the next-nearest-neighbor (frus- 
trating) exchange interaction J2 and should became rele- 
vant at some critical value J2C, past which the spin chain 
becomes dimerized and has a spin gap. Thus, the bare 
value of the coupling constant 52 should change sign at 
J2C, and close to this critical point it should have the sim- 
ple form g2 ^ a{J2 — J2c)- This will be the initial value 
for our renormalization group (RG) analysis. Instead, 
the backscattering interaction of the electronic chain is 
marginally irrelevant (for repulsive microscopic interac- 
tions) and hence the bare value of gi is negative. Thus, 



although the Hamiltonian for the spin sectors of the two 
subsystems have the same form, they are not equivalent. 

In Eq. (13.81) the Kondo term has been split into the 
53 term for the coupling of currents with the different 
chirality, and the 174 term for the coupling of currents 
of the same chirality. Their initial (bare) values are 
530 — 340 ~ uJk- Finally gi corresponds to a possi- 
ble backscattering term in the electronic chain induced 
by electronic correlations. For repulsive interactions it 
has a negative bare value, 510 < 0. 

The (chiral) spin currents of the electronic system and 
of the spin chain generate, separately, an SU{2)i Kac- 
Moody algebra. Consequently, the chiral spin currents 
•faR,L (with a — 1,2) have the operator product expan- 
sion (OPE)^ 



JUZa)jl{wt) 



SabS"^ 



8n{za - WbY 

8Tr{Za - Wh)2 



E 

7 

E 



2TT{Za - Wb) 
2Tr{Za - Wb) 



(3.9) 



where a, 13 = x,y,z; a, b — 1,2 and Za = ix + VaT, where 
r = ji is the imaginary time. 

From Eq. p.9p we can obtain the one-loop RG equa- 
tions using the procedure outlined in the Appendix A of 
Ref. 24 (or the general approach described in Refs. 22 
andl25l.) It easy to see that 54-term does not contribute 
to the lowest order RG equations. The remaining equa- 
tions are thus decoupled and read as follows 



dgi 

dl 



2'KVi 



dg2 

dl 



2nv2' 



dga 

dl 



9l 



■k{vi +V2) 
(3.10) 
Thus, starting from a negative value gi (slowly) ap- 
proaches zero, and will be neglected in the sequel. For 
next-nearest-neighbor (frustrating) exchange coupling J2 
smaller than the critical value 320 < Oj and 32 is also 
(marginally) irrelevant. In this regime the flow is con- 
troled by the marginally relevant Kondo coupling, and 
the system flows to the same fixed point of the unfrus- 
trated Kondo-Heisenberg chain. In this regime the sys- 
tem is in a PDWphase of the unfrustrated chain, as 
discussed in Ref. [3| (using results from Refs. [l|, y, and 

io 

We now follow the results of White and Affleck?- to 
analyze the situation for J2 > J2c- In this regime the bare 
value of (72 is positive, (720 > Oj and the backscattering 
interaction of the spin chain is marginally relevant. We 
thus find two stable phases depending on which coupling 
52 or 53 reaches first its strong coupling limit under the 
RG flow of Eg. p.lOp . If 52 wins, the fixed point describes 
a dimerized spin chain (with a spin gap) and a decoupled 
electronic system. In this limit the Kondo coupling is 
irrelevant. Conversely, if 53 reaches the strong coupling 
limit first, the system is again in the PDW phase in which 
dimerization is suppressed. 



We can see how these phases arise using abelian 
bosonization of the SU {2)i Kac-Moody algebra,?^. In this 
case the chiral spin currents can be represented in terms 
of the chiral bosonic fields ipsa (with s — R,L and a = 1, 2 
representing each Heisenberg chain) 



^^''^ 2TTa^ 



=piV87r'^ 



F — 

■-'sa ~ 



1 



dx^sa (3.11) 



These chiral currents have (as they should) scaling di- 
mension 1. In this representation the spin chain backscat- 
tering interaction, the 52 term in Eq. p.8|) . becomes 



ff2J 



R2 



'L2 



92 



(27ra)2 



92 



cos(v87ri^2) + -^dx'iyR2dx(pL2 



2n 



(3.12) 



where ipa == ^fiRa + fLa, and 9a = (fiLa - fRa is its dual 
field. 

In the regime in which the RG flow drives 172 to strong 
coupling, the field (^2 is pinned at the minimum of the 
cosine, i.e. at the values ip2 — {n + \)^J^. In this phase 
the order parameter e defined in Eq. (j2.5p . has a non- 
vanishing expectation value, and the system dimerizes 
(the D phase in what follows) breaking the translational 
symmetry and has a finite spin gap in its spectrum. Since 
the magnetic chain has a spin gap, its spin-spin correla- 
tion functions are short ranged. However, White and 
Affleck^ found that in the dimerized phase the spins of 
the magnetic chain have short-ranged incommensurate 
antiferromagnetic order. On the other hand, the elec- 
tronic chain in this phase remains decoupled with gapless 
charge and magnetic excitations. 

In the opposite case, in which 53 flows to strong cou- 
pling, the system flows to the pair-density-wave fixed 
point of the unfrustrated Kondo-Heisenberg chain. Let 
us summarize, for completeness, how this happens. The 
inter-chain backscattering operators have the bosonized 
expression 

Jfl.! • Jl2 + JlI • J_R2 = 



(27ra)2 



cos[V87r((p_Ri -I- ipL2)\ + cos[V87r((pLi + 'fR2)\ ) + 

dx^Rldx(pL2 + dxipLldxipR2 



1 

2^ 



(3.13) 



The first term in Eq. (|3.13|) can be written in the form 

—^1 cos[2V^((^+ + 0-)] + cos[2V^(¥>+ -0-)]\ 

(3.14) 

where (f± = ^ ((^2 ± fi) and e± = ^ (6*2 ± 6*1). 

The operator shown in Eq. p.l4p is the marginally rel- 
evant interaction found in Ref. [5] . Its presence in the ef- 
fective action drives the system to a fixed point in which 
the fields 2^/TTip^ and 2^/tt6^ are pinned at the values 



2mT and (2n + l)7r, or (2n + 1)tt and 2mr. The corre- 
sponding phase was analyzed in Refs. p] and [7] where 
it was shown that it has gap for ah spin excitations. In 
this phase there stih is a decoupled gapless charge sector. 
A remarkable feature of this phase is that the only or- 
der parameters that exhibit quasi-long-range order (i.e. 
have power-law correlations) are composite operators 
made from observables of the electronic chain and the 
spin chain that, separately, have short-range correlations 
since their correlation functions fall-off exponentially fast 
with distance. In particular this phase is characterized 
by the pair-density wave order parameter. 



^PDW 



^TS 



•AT, 



(3.15) 



where Ats = iY.a,fi^aR,i{x)[cTay)a^p^aL,i is the spin 
triplet pairing operator of the electronic system ( "chain 
1" ), and N2 is the Neel order parameter of the spin chain 
("chain 2"). Thus, the PDW order parameter is a four 
fcrmion operator. The PDW operator of Eq. p.l5|) has 
the bosonized form 



A 



PDW — 



2(7ra)2 
cos(-\/47riy9_)] 



[2 cos(v47r6'_) + cos(v47r(^+) 

(3.16) 



where we have omitted an oscillatory factor with wave 
vector QpDW — ^r. From this expression it is appar- 
ent that in the PDW phase the operators in brackets 
in Eg. p. 161) have finite (and non- vanishing) expectation 
values and, as a result, the PDW order parameter has 
power-law correlations of the form 



^PDW 



{x)i^l^^m 



1 



j^flPD'W 



(3.17) 



The exponent takes the value r}PDW = 1 if the electronic 
chain is non-interacting, while for repulsive interactions 
(with Kc > 1) it increases to the value ryp^n^ = K^ > 1. 
Berg et alM^ showed that an ordered PDW state with 
ordering wave vector Q always has a subleading uniform 
superconducting order but with charge 4e instead of 2e as 
in a conventional superconductor. In Ref. ;7] showed that 
the spin-gap phase of the Kondo-Heisenberg chain has 
PDW quasi-long-range order as well as charge 4e quasi- 
long-range order (albeit with a larger critical exponent). 
For a PDW with ordering wave vector Q = n (as in 
the present case), the charge-4e superconducting order 
parameter A^e is simply the square of the PDW order 
parameter. Hence we can make the identification 



const. X e 
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(3.18) 



This order parameters has scaling dimension 2Kc, and its 
correlation function falls-off with an exponent 7]^^ — "iKc- 
The RG fiows of Eq. (|3.10l) determine the structure of 
the phase diagram in this weak coupling regime. The RG 
flows are marginally unstable on both coupling constants 
and hence the effective couplings run to their strong cou- 
pling regime. As usual, the velocities are only affected by 



irrelevant operators and acquire at most a finite renor- 
malization. Hence the velocities do not affect the RG flow 
of the dimensionless coupling constants 32 and 33 since 
gi is a marginally irrelevant coupling. The dimensionless 
coupling constants have initial values 320 ~ (^2 — J2c)/t 
and gs -- Jn/t. 

At the level of the one-loop RG of Eq. (|3.10p . which 
is accurate for 32 and 173 small, the location of the 
phase boundary (the separatrix of the RG flow) be- 
tween PDW and the dimerized phases is the straight line 

(73 = Y' ( 1 + ^ ) or, equivalently in terms of the macro- 
scopic parameters, Jk = 5(^2 — J2c) (l + xj- Above 
this phase boundary the system is in the PDW phase, 
and below this phase boundary it is in the dimerized 
phase. Along the phase boundary both coupling con- 
stants §2 and (?3 flow to strong coupling at the same rate. 
Thus, the quantum phase transition between the PDW 
phase and the dimerized phase is likely to be flrst order. 
We conclude that quantum frustration of the spin 
chain leads to irrelevant and essentially unobservable ef- 
fects in the PDW phase. On the other hand quantum 
frustration leads to a dimerized phase and the Kondo 
coupling has essentially irrelevant effects in this phase. 
This is so because in this phase there is a spin gap in 
the (frustrated) spin chain producting an effective de- 
coupling from the electronic chain. It is well known (and 
easy to see) that if the spin chain is treated classically, 
its ground state is an incommensurate spiral. The nu- 
merical (DMRG) calculations of White and Affleck show 
that there is short-range incommensurate order inside the 
dimerized phase. 



IV. THE J2 > Ji REGIME 

We now turn to the limit in which the next-nearest- 
neighbor exchange interaction J2 is strongest than Ji . In 
this regime the treatment used in the preceding section 
is not adequate. Instead, it is convenient to think the 
magnetic chain as a zig-zag two-leg ladder, i.e. two spin 
chains with nearest neighbor exchange J2, and weakly 
coupled by a inter-chain (zig-zag) exchange interaction 
Ji. In this picture the frustrated spin chain is a trian- 
gular ("trestle") ladder. In the extreme case Ji = the 
system we are considering becomes two Heisenberg an- 
tiferromagnetic chains coupled with an electronic chain 
by the Kondo exchange interaction. In the limit in which 
the Kondo coupling is also weak, the starting point of our 
treatment consists of three chains (two magnetic chains 
and an interacting electronic chain) weakly coupled by 
Jk and Ji . This regime turns out to be quite rich. 

The zig-zag antiferromagnetic ladder in the limit J2 3> 
Ji was considered by White and Affleck^ who treated 
the system as two weakly coupled Heisenberg chains. In 
this limit, the effective low energy theory of this zig-zag 
ladder includes a marginally relevant backscattering cou- 
pling between the SU(2) spin currents of the two chains. 



However, due to the special symmetry of this ladder the 
(otherwise relevant) inter-chain coupling of their Neel or- 
der parameters is absent. White and Affleck also showed 
that the marginally relevant inter-chain backscattering 
interaction causes the RG to flow to a strong coupling 
fixed point with finite (albeit exponentially small) dimer- 
ization and an also exponentially small spin gap. These 
results were confirmed using a DMRG. Furthermore in 
this regime the zig-zag ladder has short-range incommen- 
surate spiral order. Thus, for the problem of interest 
here, we conclude that if the Kondo exchange interac- 
tion is also very small, the low energy physics is that of a 
dimerized zig-zag chain essentially decoupled from a one- 
dimensional Luttinger liquid with gapless and decoupled 
charge and spin excitations. This is the same phase that 
we encountered in the preceding section. 

Another limit of interest is the case in which Ji — > 
while holding Jk finite (and small) , in which the system 
reduces to two antiferromagnetic chains and a conducting 
fermionic chain coupled by a Kondo exchange interaction 
Jk- Up to some simple redefinitions, in this limit the 
problem we are interested in is related to the problem 
of over-screened (or mult i- channel) Kondo-Heisenberg 
chains considered by Azaria and Lecheminant,^^ (who 
also considered the frustrated three-leg Hciscnbcrg lad- 
der.) We will see below that the spin sector of the prob- 
lem we are interested in is equivalent, in the limit Ji — ^ 0, 
to the spin sector of an over-screened Kondo-Heisenberg 
model. 

The main result of the work by Azaria and Lechem- 
inant is that the Kondo exchange interaction, which is 
marginally relevant at J2 = 0, drives the spin sector to 
a finite infrared- stable fixed point with non-trivial prop- 
erties. This non-trivial fixed point is a one-dimensional 
analog of the multi-channel Kondo problem. At this fi- 
nite fixed point the spin sector decouples into two "chi- 
rally stabilized" spin liquids, first discussed by Andrei, 
Douglas and Jerez3i We will use the exact solution of 
Azaria and Lecheminant for the two-channel case, which 
is relevant to our problem, and discuss the effects of the 
coupling Ji at this non-trivial fixed point. 

To see how this works let us consider the effective 
(bosonizcd) low energy field theory of the frustrated 
Kondo-Heisenberg chain in the limit J2 ^ Ji and 
•h ^ Jk- In the low energy limit the frustrated Kondo- 
heisenberg chain system exhibits spin-charge separation. 
The total effective low energy Hamiltonian density Ti, — 
He -I- "Hs is a sum of the Hamiltonian He for the charge 
sector and Hg for the spin sector. 

The charge sector is described, as in the preceding sec- 
tion, by the charge Bose field (j)e and its dual field 9e, 
with a fixed charge Luttinger parameter Ke and a ve- 
locity Ve- The effective Hamiltonian for the charge sec- 
tor. He, describes the charge degrees of freedom of the 
fermionic system and, hence, is the same as the one given 
in Eg. (1X71). 

In this regime, the spin sector is described by the chiral 
spin currents of the two weakly coupled legs of the zig-zag 



ladder and the chiral spin currents of the ID electronic 
system. The Hamiltonian density Hs of the spin sector 
is given by 



Hs —Hq + Hint 

no — — - — 



'R.A 



"fl,l- 



Jl,i J 



L,l- 



'■ Jit',2 • J_R,2 : 

2nv2 



Jl,2 • Jl,2: 

Jfl,3 • Jfl,3 : + '- Jl,3 ■ Jl,3 : 



(4.1) 



Hint —gi 



(Jfl,,l + Jfl,2) ■ Jl,3 + (Jl,1 + Jl,2) • J_R,3 



+52 (J_R,1 • Jl,2 + J_R,2 • Jl,i) 

+9 (Jfl.l • Jl,1 + J_R,2 • Jl,2) + 9 Ji?,3 • Jl,3 

(4.2) 

Here 1 and 2 label the SU(2) spin currents the two mag- 
netic chains and 3 labels the SU(2) spin current of the 
electronic chain. The bare values of the spin velocities, 
vi and V2, and of the effective coupling constants gi 
and (72, are respectively given by vi ~ aJ2, V2 — at, 
gi ~ aJK, 92 — aJi. The last line in Ea. (|4.2p represents 
the marginally irrelevant backscattering intra-interaction 
of the spin currents of the two spin chains. The effective 
coupling constants g and g' are given by: g ~ — Ji < 
for the spin chain, and g' w — 2[//t < for the ID elec- 
tronic system (where U and t are the on-site Hubbard 
interaction and t is the nearest-neighbor hopping ampli- 
tude), respectively. As in Section iHll in Eq. (|4.2|) we 
have neglected the redundant term that couples the spin 
currents with the same chirality, which leads to a finite 
renormalization of the velocities. 

Also as in Section [TTTl here too we have not included in 
Eq. (|4.2|) a coupling between the spin-density wave or- 
der parameter of the ID electronic system and the Neel 
order parameters of the weakly coupled spin chains since 
their associated ordering wave vectors are different and 
hence the coupling is not allowed in a system with trans- 
lation invariance. Similarly, we have not included the 
(relevant) coupling between the Neel order parameters 
of the two weakly coupled spin chains since it is forbid- 
den by the symmetries of the zig-zag chain. The absence 
of these (potentially most strongly relevant) interactions 
has important consequences for the stability analysis of 
the fixed points and for the structure of the phase dia- 
gram. 

The effective field theory for the spin sector described 
by Eqs. ((4T|l and (|42|) has three SU(2)i spin cur- 
rents (two SU(2)s for each spin chain and one SU(2) 
for the electronic system). Thus it is a perturbed 
SU(2)ixSU(2)ixSU(2)i conformal field theory. In ad- 
dition the system has a U{1) charge sector which (in the 
low energy limit) is decoupled from the spin sector. In 
the absence of the perturbations of Eq. (|4.2p the spin sec- 
tor has central charge c = 3. 
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A. Jk <^ Ji <^ J2' Dimerized Phase 

The effect of the perturbations is to open up gaps thus 
driving the system to a fixed point with a smaller central 
charge. The RG equations for the effective field theory 
of Eq. (|4.2p have the same decoupled form as before, i.e. 



dgi 

dl 



9l 



■n{vi +V2)'' 



dg2 
dl 



2'KV^ 



(4.3) 



The RG flows are decoupled, as in the case of section Hill 
Which gap opens up depends on the relative strengths of 
the perturbations. 

For Jk <C Ji, 52 ^ ffi, the RG flows to the strong 
coupling fixed-point (of .92) of the dimerized zig-zag chain 
which is now gapped. Since the perturbation with cou- 
pling constant 32 is marginally relevant the gap is expo- 
nentially small (up to logarithmic corrections due to the 
marginally irrelevant coupling g)^ In this phase the ID 
electronic system and the zig-zag chain are decoupled at 
low energies. This is the same dimerized phase we found 
in Section Hill with the only difference that when J 2 ^ Ji 
the zig-zag chain has short-range incommensurate spiral 
order. ^ Thus, in this phase we have a gapless charge mode 
and a gapless spin mode, both belonging to the decoupled 
electronic system. 

On the other hand, for Jk ^ Ji the RG flows to differ- 
ent strong coupling regime that will be described below. 
From the results of Section IIIII we know that one possi- 
bility is that this phase may also be the PDW phase that 
arises for Jk '> Ji (and J2 <C Ji). The other possibil- 
ity is that the strong frustration regime is in a different 
phase altogether. 

The phase boundary that separates the dimerized 
phase from the non-trivial phase is the separatrix of 
the RG flow of Eq. (|4.3l) . This fiow has a separatrix at 

1 + ^] or, equivalently, "^ = 5^ (l + :^)- 



.91 



£2 
2 



B. The Ji < Jk < J2 regime 

The nature of the phase in the Ji <C Jk ^ J2 regime 
cannot be accessed by the perturbative RG of Ea. (l4.3p 
and must be determined non-perturbatively. This regime 
is governed by a finite, infrared stable, fixed point that 
was first discussed by Azaria and Lecheminan t^^'^^ (who 
were interested in a somewhat different problem) who 
constructed a non-trivial fixed point of a type first dis- 
cussed by Andrei, Douglas, and Jerez^*^ using a non- 
perturbative approach. We will consider first the case 
Ji — and use the results of Refs. [15] and 18] to con- 
struct the non-trivial fixed point at a finite value of the 
Kondo coupling constant. Once this fixed point is identi- 
fied we will assess the role of the perturbation with cou- 
pling constant Ji at the non-trivial fixed point in order 
to determine the actual nature of the phase. In particu- 
lar if the fixed point is unstable even for small Ji we will 
see that the resulting phase has the same properties as 



the pair-density-wave phase of the unfrustrated Kondo- 
Heisenberg chain. However if the non-trivial fixed point 
is stable then there is a new phase (for weak enough 
Ji) which is distinct from the PDW phase and from the 
dimerized phase. 



1. Ji = 0: Non-tnvial Finite Fixed Point 

We begin by considering first the case Ji — and 
summarize the construction of the nontrivial finite fixed 
point by Azaria and Lecheminant.— So, we will consider 
the spin sector of our system with Hamiltonian T-Lg, see 
Eas. (|4.ip and (|4.2p . and set for now g2 = g = g' = 0. 
The Hamiltonian density now reads 



ns = 



2'KVl 



y ^ ( : J(,fl • Ji,_R : + : Ji,L • Ji 



j=l,2 

+ — T^l : Js.R- J 



o ( ■ •^3,R ■ Ja.i?, : + ■ Js.L • Js.L ■ 

+91 ( Y. J^^^ ■ J3,L + Yl -^^.^ ■ *^3,fi) (4-4) 

'i=l,2 1=1,2 

The free part of this Hamiltonian has an SU(2)iX 
SU(2)iX SU(2)i symmetry which is partially broken by 
the interaction term. Since the interaction retains the 
symmetry of the exchange of the spin sectors of the two 
spin chains (here denoted by 1 and 2), it is natural to 
rewrite this system in terms of the total chiral right- and 
left-moving chiral currents current Ir^l — Ji.i? + J2,fl- 
The chiral currents Ir^l are the generators of two chiral 
SU(2))2 Kac-Moody algebras (one for each chirality). 

We will first rewrite the two decoupled chains in terms 
of a theory of four Majorana fermionsi22i2^ This is pos- 
sible since a conformal field theory with current algebra 
SU(2)iX SU(2)i is equivalent to a conformal field theory 
of the Wess-Zumino-Witten model with current algebra 
SO (4)1. This model has a well known representation in 
terms of four free Majorana fermions^S 



27ri;i 



-E(^J^«-J 



■ Ji,fi 



Jj r, • J 



L ■ •li,L 



2 



4=0,1,2,3 

where £,l{^) (with i = 0,1,2,3) denotes the four 
species of Majorana fermions with both chiralities (s — 
R, L), and satisfy canonical anti-commutations relations, 
{e(a;),e} = S.jSix-y) and {?l,(x),ei(y)} = 0. 

In terms of the Majorana fields, the chiral SU(2)2 cur- 
rents, Ig (with a = 1, 2, 3 and chirality s — R,L), which 
are three of the six chiral currents of S0(4)i (with both 
chiralities), are given by 



i: = JI. + J2 



^abc /-b /-c 



CsCs 



(4.6) 



The remaining three chiral currents of S0(4)i are (with 
s = R,L) 



K[ 



=^ces 



(4.7) 



where K^^^l = Ji,fl,L — J2,fl,L- This decomposition 
allows us to rewrite the two decoupled spin chains in 
terms of the SU(2)2 chiral currents. It is equivalent to 
the identificatioitS^i^i^ SU(2)iX SU(2)i ~ SU(2)2 x Z2, 
where the Z2 factor will be represented by a critical Ising 
model, i.e. a massless Majorana fermion which cannot 
be written in the Sugawara form (as a quadratic form in 
chiral currents). 

Following Ref. [l5| we now rewrite the Hamiltonian 
density of Eq. (|4.4p in the form 

ns = -y (efl*5,C° - eLid^^L) + Hi + H2 (4.8) 

where 

rii ——— -.Ir-Ih: H — - — : J^^l ■ J3.L ■ +3iI_R ■ Ja.L 

TTUl 27ri;2 , , , T T 

rL2 ——z- ■ i-L ■ i-L ■ H — :; — : ^3,^; • Js.fl '■ +9ii-L ■ •13, b. 

(4.9) 

Obviously, "Hi and 'H2 commute with each other, 
[Hi,H2] — 0. The Hamiltonians "Hi and H2 involve 
the chiral SU(2)2 currents 1b,l and the chiral SU(2)i 
spin currents J3,fl,L of the electronic system. In addi- 
tion, these two Hamiltonians also commute with the free 
Majorana Hamiltonian shown in Ea. (l4.8p . which repre- 
sents a Z2 (Ising) sector. The Hamiltonians "Hi and 7^2 
have the form of the chirally-stabilized spin liquid theory 
of Andrei, Douglas and Jerea^S who showed that each de- 
scribes a theory controlled by a finite infrared-stable fixed 
point. Although the fixed point of the theory of Hamil- 
tonian Hi is not chirally-symmetric, the sum "Hi +^2 is 
chirally-symmetric. 

The most direct way to identify the finite non-trivial 
fixed point is to use abelian bosonization approach of 
Azaria and Lecheminant^'^ which solves the problem by 
formally breaking the SU(2) symmetry down to a U(l) 
subgroup. In order to do this we will formally break the 
SU(2) symmetry of the interaction terms of the Hamil- 
tonians Til and T-L2 which will now read 

"Hi^nt ^giwInJ^L + ^ {^R'h + ^R-hh) 

H2,i„t ^gnllJiR + ^ {itJ^R + IZ4r) (4.10) 



with gin = gi± = gi in the SU(2)-invariant case. The 
resulting Hamiltonian is solvable at a particular point, 
known as a Toulouse point, upon a simple unitary trans- 
formation. This approach is well known from the the- 
ory of the Kondo problem, and was used with great suc- 
cess by Emery and Kivelson in the two-channel Kondo 
problem. -^^ The SU(2) invariant Hamiltonian of Eas. (|4.8p 
and (|4.9p (solved by Andrei et al. using Bethe ansatz 



methods^) and its U(l)-invariant version are also related 
by an irrelevant operator at the Toulouse point. 

A summary of the Azaria-Lecheminant solution at the 
Toulouse point is presented in Appendix[Xl The upshot of 
this solution is that the low energy degrees of freedom at 
the Toulouse point are a massless boson, that we denote 
by $1, and two massless Majorana fernuons, denoted by 
^° and x^ respectively. Since the Bose field $1 is at the 
SU(2)i radius, the Luttinger liquid Hamiltonian Hll is 
equivalent to a theory of an SU(2)i WZW model with 
right- and left-moving chiral currents ^r and J'l 



-7± _ „^iVSTr'S>l.R 



rr± _ _|_ ±i\/87r$i 

^^ " 2^a 



^^.= 



Jl = 



1 



/2^ 
1 

7^ 






(4.11) 



In the notation of Appendix [Aj the effective Hamiltonian 
at the Toulouse point is 

Htp - ^{Jl + Jl) +nF[e]+'HF[x'] (4.12) 

where "HFiC^] and HFix^] are the Hamiltonians for the 
free Majorana Fermi fields ^° and x^, given in Eq. (jA15[ ) 
and Eq. (|A16[) of Appendix [X] The effective low en- 
ergy Hamiltonian of Ea. (|4.12p is the chiral fixed point 
of Azaria and Lecheminant. It describes a conformally- 
invariant system with total central charge c = 2 (in the 
spin sector alone), down by one unit with respect to the 
decoupled system. 



2. Gapless Fractionalized Spin Fluid 

To assess the significance of this fixed point to the frus- 
trated Kondo-Heisenberg chain we need to consider the 
stability of this fixed point. The results of Azaria and Le- 
cheminant will be useful but require some changes given 
the differences in the two problems. In particular they 
showed that the corrections to the Toulouse point Hamil- 
tonian that correct for the artificial breaking of the SU(2) 
symmetry amounts to an irrelevant perturbation. 

We will now see that for Ji weak enough there is a 
novel gapless fractionalized fluid phase characterized, in 
the spin sector, by the stable fixed point of Azaria and 
Lecheminant. In this phase we will find that in the charge 
sector there is coexistence between the PDW order with 
the singlet superconductor (and other more conventional 
orders) . However both in the charge and spin sectors the 
leading operators, which have the smallest scaling dimen- 
sion and hence the more strongly divergent susceptibil- 
ities at low temperatures, are unconventional. In this 
sense this phase is characterized for having intertwined 
orders in a fractionalized fiuid state. 

As we stressed before in the problem at hand we do not 
have a coupling between the Neel order parameters of the 
spin chains with each other and with the ID electronic 
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system. We thus need to only consider the role of the 
backscattering interactions between the two spin chains 
(with coupling constant 32 in Eq (|4.2|) ). and the backscat- 
tering interactions in each spin chain (with coupling con- 
stant g in Ea. (|4.2p ') and in the ID electronic system (with 
coupling constant g' in Eq. ()4.2|) ). We wiU further assume 
that the coupling constants g and g' are approximately 
equal to each other. The results of Ref. l15|] imply that 
in the low-energy limit both the inter-chain spin-current 
backscattering interactions and the intra-chain backscat- 
tering spin-current interaction (including the ID elec- 
tronic system) map onto the operator jTfl, • JTl (although 
with different effective coupling constants) where ^r^l 
are the chiral SU(2)i currents defined in Eq. ()4.11|) . 

Given these considerations we now can write the full 
Hamiltonian for the charge and spin sectors H = He + 
Hs, where He is the Hamiltonian for the charge boson 
(j)c presented in Eq. (|3.7|) . and Hs is the Hamiltonian for 
the spin sector of Eq. (|4.2|) at the chiral fixed point (the 
Toulouse point) 



27I'-Ul I 2 



3 
2 



R -r- J L + QcsJr ■ J L 



{xrQxXr - x\dr.x\) 



(4.13) 



where g^s = 3^ -f 2*72 is the effective coupling constant 
for backscattering interactions of the SU(2)i currents of 
Eg. (14. Ill) . This interaction with coupling constant gcS^ 
which in abelian bosonization leads to an operator of 
the form cos(-\/87ri'i) in the effective Hamiltonian (see 
Eg. (14. Ill) '), is a marginal operator. It is a relevant per- 
turbation for (7eff > and an irrelevant perturbation for 
f/off < 0. Since g2 ^ Ji > ^ and g ^ —Ji < 0, the rel- 
evance or irrelevance of this interaction depends on the 
relative strengths of these couplings. In the limit that 
we are considering here gcs < until the interaction Ji 
becomes large enough. 

We then conclude that, at least for weak enough in- 
terchain interaction Ji, the chiral fixed point is pertur- 
batively stable provided g^s < 0. Hence in the regime 
Ji ^ Jk ^ J2 there is a new stable phase which is char- 
acterized by the finite fixed point, which is described by 
the Toulouse Hamiltonian of Eq. (I4.13P (up to irrelevant 
operators). We will now look at the behavior of the corre- 
lators of the frustrated Kondo-Heisenberg chain. To this 
end we need to find the form of the physical observables 
of the Kondo-Heisenberg chain in this description. 

Azaria and Lecheminant gave an identification-'^^ at 
the Toulouse point for the Neel order parameters Ni 
and N2, as well as the SDW order parameter of the 
ID electric chain N3, and of the spin currents Ji, J2 
and J3. Here we will adapt these methods to identify at 
the Toulouse point, in terms of the fields of Eg. (14. 131) . 
the observables of the ID electronic system Ass{x) and 



^Ts{x), the singlet and triplet superconductor oper- 
ators, and its CDW and SDW ("Neel) order parame- 
ters pcDw{x) and Nsdw{x) (both with ordering wave- 
vectors 2kp). We will also need the Neel order parame- 
ters Ni and N2 of each magnetic chain (both with order- 
ing wave- vector tt), and in the composite order parame- 
ters Ats ■ (-^1 ±-^2), which describe a pair-density-wave 
order parameter A^^^ also with wave-vector tt. As in 
the case of the unfrustrated Kondo-Heisenberg chaiu;^ 
the PDW order parameters involve the degrees of free- 
dom of the ID electronic system and of the antiferromag- 
netic chains. 

The spin currents J3 and the SDW (Neel) order pa- 
rameter N3 of the ID electronic chain at the Toulouse 
point are identified (to leading order) with the following 
operator 



,15 



J. 



3,B. 






JL 


7± ^ 


-Jt 


^Ri 




■^Ji 


*™xl?xiA/'^ 


^3"^ 


- *7ra XrxI^^ 

(4.14) 



where the fields N^ and M^ are given by 



A/"^ sin(V2^$i) 

AA± _±e±^V2^ei 
■na 

8 ^ — cos(V27r#i) 



(4.15) 



are the components of the (spin-1/2) primary field g 
(with scaling dimension 1/2) of the SU(2)i WZW con- 
formal field theory 



£ + iAf -a 



(4.16) 



where cr are the three Pauli matrices. 

The spin currents Ji and J2 of the spin chains, and 
the associated Neel order parameters Ni and N2 have 
the following (leading order) operator identification at 
the Toulouse pointiS. 



jfR-jRX^+^a^xieR 

J2R-ji{}-^0.^x\eR 

-fZ JZ rrZ 



(4.17) 



where a is the short-distance cuttoff (the lattice spacing). 
Similar expressions hold for the left moving components 
(with the R and L labels exchanged). 

The Neel order parameters Ni and N2 have the fol- 
lowing identificationai^ 

Af -A/'^(/i5/io + 0-5(70) 

A|-7V^(/l5/10-cr5CTo) 

A^f*" --A/'='= (/i5Aio - ctso-q) 



^2 ^^f^{fJ,5^J■0 + cr5Cro) 



(4.18) 
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Here ao and us, and /xq and fj,^ are respectively the or- 
der and disorder twist field operators associated with the 
critical Ising models with Majorana fields ^° and x^ (see 
Appendix [XJ . 

We will now turn to the identification at the Toulouse 
point of the order parameters involving the charge sec- 
tor of the frustrated Kondo-Heisenberg chain. They are 
the CDW order parameter pcDW of the ID electronic 
system, 



1 



PCDW xWl e-^^^."^. cos(V2^|.i) (4.19) 

na 

the singlet and triplet superconductor order parameters 
Ass and Ats of the ID electronic system, 



A 
Ats 



1 



SS XrxI e-*^2.e. cos(V2^$i) 

na 



-xWl e-'^'^ M 



(4.20) 
(4.21) 



and the PDW order parameters the the frustrated 
Kondo-Heisenberg chain A^^^ = A • (Ni ± Ng), 



^- ^ M5M0 e 



(4.22) 
(4.23) 



where we used that the Majorana mass term ix\ix\ ^ 
ei where ei is the energy operator of the Ising model 
(the scaling dimension 1 relevant "thermal" operator at 
the critical point) and the operator-product expansion 
(OPE) 0-5 El ^ (75 (likewise the disorder operator satisfies 
Ms^i ~ Ms)- Similarly the Majorana bilinear 1^7??° is 
identified with the energy density operator e° of another 
critical Ising model and satisfies the same OPEs with its 
order and disorder operators, (Tq and /iq- 

Finally, the charge-4e uniform superconducting or- 
der parameter A4e (already discussed in Section IIIIl 
Eg. p. 181) ) which is the square of the PDW order pa- 
rameters (given in Eq. (|4.22p and Eq. (|4.23p ). is in turn 
identified with A4e ^ exp(i2-\/27r6'c). This operator has 
scaling dimension 2Kc- Notice that here we used the fact 
that to leading order the Ising twist fields a and p, to lead- 
ing order fuse into the identity operator. It is easy to see 
that the next-to-leading term involves the energy density 
operator e of the Ising model. Since the energy density 
operator e has scaling dimension 1, the leading correc- 
tion to the PDW order parameter has scaling dimension 
4i^, + l. 

We can now characterize the phase controlled by the 
non-trivial (finite) chiral fixed point represented by the 
effective Hamiltonian of Eq. (|4.13p in the regime gcff < 0. 
At this fixed point (and in the low energy regime of the 
entire phase which it controls) the four fields (f)c, ^1, <^° 
and x^ are massless and hence critical. The only effect 
of the marginally irrelevant perturbation J' ^ ■ jTi is to 
induce logarithmic corrections to the correlators (except 
those that involve conserved currents). 



a. Order parameters of the spin sector : The scaling 
dimensions of the spin currents 3k,R,L (with k = 1,2, 3) 
is 1. The scaling dimension of the SDW order parameter 
N3 (given in Eq. (|4.14p ) is | since this operator is the 
product of the Ising energy density operator ei (which 
has scaling dimension 1) and of the SU(2)i primary field 
A/" (which has scaling dimension i). Instead, the Neel 
order parameters Ni and N2 (given in Eq. (|4.18|) of the 
spin chains is 3/4 since each disorder (and order) operator 
has dimension i In contrast, at the decoupled fixed 
point of Eq. (|4.2p with gi — 0, the scaling dimensions of 
all three Neel order parameters is ^. Thus, the Kondo 
exchange coupling has caused the scaling dimensions of 
the three Neel order parameters to increase considerably 
and the critical exponents of their correlation functions 
are now r/i = 772 = 3/2 and 773 = rjsDW = 3. 

We should also consider the spin singlet composite op- 
erators 9T± = N • (Ni ± N2). Both operators have order- 
ing wave vectors Q± — ■n~2kp. Since these operators are 
spin singlets, they can be interpreted as composite CDW 
order parameters. A similar operator also exists in the 
weakly frustrated case^i In addition, the operator VI- is 
odd under the exchange of the two magnetic chains (as is 
also the odd PDW order parameter /\P^DWy ^t is easy to 
see that at the Toulouse point the operators 91± are iden- 
tified (to leading order) with the operators 01+ ~ pzPo 
and Ol_ ^ CTsCToi respectively. Hence, both operators *n± 
have scaling dimension ^ and hence their critical expo- 
nent is rj<ji_^ = 2- Therefore in this phase, the dominant 
order in the spin sector is given by the composite order 
parameter 0T± , followed by the Neel orders of the mag- 
netic chains and by the SDW order parameter of the ID 
electronic system. 

b. Order parameters of the charge sector : We can 
also read-off that the scaling dimensions of operators that 
involve the charge sector. The scaling dimensions of the 

CDW order parameter is 5 ( 3 -I- -j^ ) at the Toulouse 
point, while the scaling dimension of the singlet super- 
conducting order parameter is \{Kc + 3) (as it is for 
the triplet superconductor), and the scaling dimensions 
of the PDW orders is ^ i^Kc + \) and for the charge-4e 
superconductor is 2Kf.. Their critical exponents are, re- 
spectively, given by 'qss = Vts = ^c+3, Vcdw =3+;^, 
Vpdw ^ Kc + ^ and 7/4e = AKc- 

From this analysis we conclude that in this phase all 
of the orders are present, with the composite CDW or- 
der parameters *yi± and the PDW superconducting order 
parameter Apow having the correlation function that 
decay more slowly with distance, and hence are the dom- 
inant orders in this phase. In this sense this phase is 
analogous to the PDW phase of the weakly frustrated 
regime. However, in this phase these unconventional or- 
ders also coexist with the more conventional Neel orders 
of the chains and in the electronic system, as well as the 
singlet (and triplet) superconducting orders, the CDW 
order and the charge 4e superconducting order. There- 
fore this phase has a gapless charge sector (with central 
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charge c = 1) and a partially gapped spin sector (with 
total central charge c = 2). This fractionalized phase has 
a large number of orders that are comparably strong and 
hence are intertwined (instead of competing with each 
other). This is a feature that is shared with the PDW 
phase of the weakly frustrated regime. In contrast, in the 
dimerized phase the frustrated Heisenberg chain and the 
ID electronic system are simply decoupled at low ener- 
gies, and the allowed orders are more conventional. 



3. Quantum phase transition to the PDW phase 

We will now consider the case in which g^fi > in 
which case the operator J'jj^ ■ J'l becomes marginally 
relevant. This is a Kosterlitz-Thouless type transition. 
In this phase the field $i becomes pinned at the val- 
ues (n 4- ^) -y/^, and its fluctuations are massive. Con- 
sequently there is a gap in most (but not all) the spin 
degrees of freedom. This leads to a further reduction of 
the central charge of the spin sector of at least from c = 2 
to c = 1, while the charge sector remains massless and 
has central charge c = 1. We will see however that the 
spin sector is fully gapped and that the resulting phase 
is equivalent to the PDW phase we discussed in the weak 
frustration regime. 

In the regime in which ^eff > 0, operators such as 
cos(v27r$i) have vanishing expectation values in this 
phase, while operator such as sin(-\/27r$i) do not. More- 
over, in this phase the dual field Qi fiuctuates wildly 
and vertex operators of the dual field have vanishing ex- 
pectation values. Therefore in this phase the field $i 
is massive and decouples from the low-energy physics. 
The effective Hamiltonian for this phase has the same 
form as the effective Hamiltonian at the Toulouse point, 
Ea. (l4.12p . except that the SU(2)i currents ,J'r^l have 
been projected out. Hence, in this phase the remaining 
massless degrees of freedom are the charge boson 0c and 
the Majorana fermions ^° and x^- We will see, however, 
that these Majorana fermions actually become massive. 

These considerations allow us to determine the behav- 
ior of the observables of interest. It is easy to see that 
all the order parameters that are not expressed as com- 
posite operators have exponentially decaying correlators 
and are not condensed. This includes, the singlet and 
triplet superconducting order parameters A55 and A^s, 
the SDW and CDW order parameters N3 and pcDW of 
the ID electronic system, and the Neel order parameters 
Ni and N2 of the spin chains (see Eqs. (I4.14|) . (|4.18l) . 
(HTTO)) . (|T^ and (|OT|) '). All these operators have ex- 
ponentially decaying correlation functions. Nevertheless 
several composite operators have power-law correlation 
functions. 

a. Order parameters of the spin sector : The only 
operators from the spin sector that have power-law cor- 
relations in this phase are the spin-singlet operators *Jt± . 
As we saw, these operators are given by products of Ising 
twist operators. Hence also in this phase these oper- 



ators have scaling dimension 1/4 and critical exponent 
T](yi^ =1/2, and their correlation functions oscillate with 
wave vector tt — 2k p. 

h. Order parameters 0/ the charge sector : The only 
operators in this phase with power-law correlations in the 
charge sector are the PDW operators of Eqs. (|4.22[) and 
(I4.23p . and the charge 4e uniform superconductor. Since, 
as we will see below, the Majorana fermions (^ and x^ 
become massive, in this phase the scaling dimensions of 
/^PDW ^^^ ^^^ g^j-g ^Yie same as in the PDW phase. 

What is the difference between this regime and the 
PDW phase? So far looks almost identical to the PDW 
phase of the weakly frustrated Kondo-Heisenberg chain, 
albeit with different critical exponents if the Majorana 
fermions ^° and x^ remain massless. Microscopically, the 
main difference between this phase and the PDW phase 
of the weakly frustrated regime is that the latter phase in 
the symmetry of exchanging the two spin chains is bro- 
ken, whereas superficially here it seems to be unbroken. 

Let us examine this question more closely. From the 
structure of the effective Hamiltonian, and of the observ- 
ables, we see that what distinguishes this possible phase 
from the PDW phase is whether the product of Ising twist 
fields crop's has (or has not) a non-vanishing expectation 
value. This is achieved by combining both Majorana 
fermions into a single Dirac fermion V' = ~m{^^ + *X^) 
and generating a Dirac mass term or, equivalently, a Ma- 
jorana mass for both Majorana fermions, i.e. the energy 
density operators of the two Ising models. In particular, 
in the phase in which {oqu^) ^ 0, the disorder opera- 
tors have vanishing expectation value, (/J.0/^5) = 0, and 
vice versa (by Kramers- Wannier duality). Consequently 
across the phase transition the symmetric PDW oper- 
ator, A;^^^ has short-range correlations, whereas the 
operator /S.^^^ becomes the PDW order parameter of 
the weakly frustrated case. 

However, for the expectation value of the product of 
twist fields to be different from zero the the two remain- 
ing Majorana fermions (^ and x^ must become massive. 
Such a mass term cannot be generated in the fractional- 
ized fluid phase. However, a term of this form is gener- 
ated from nominally irrelevant operators (whose scaling 
dimension at the Toulouse fixed point is 3 and larger) 
once the operator J" ^ ■ J" ^ becomes marginally relevant 
and the field $1 becomes massive. That this is correct 
can be seen form the identification of the spin currents, 
e.g. Ea. (|4.17p . These operators exist since at the micro- 
scopic level the symmetry of exchanging the two mag- 
netic chains is explicitly broken. However, the operators 
that break this symmetry are irrelevant at the Toulouse 
point and hence do not destabilize the gapless fraction- 
alized spin fluid phase. Irrelevant operators that break 
the symmetry of a system at a fixed point arc known as 
dangerous irrelevant operators. Such operators do not 
change the universality class of the phase transition but 
change the nature of the resulting stable phase. More- 
over, the backscattering current interactions of the spin 
currents of the two magnetic chains and of the ID elec- 
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tronic system, Eq. (j4.2[) , contain operators that break this 
symmetry but which are irrelevant at the Toulouse point 
and at the quantum phase transition out of the fraction- 
alized gapless spin fluid phase. Thus, the resulting phase 
is a PDW phase due to the existence of these dangerous 
irrelevant operators which are irrelevant at criticality but 
which lower the symmetry of the stable phase. 



V. PHASE DIAGRAM AND CONCLUSIONS 

In this paper we discussed the structure of the phase 
diagram of the frustrated Kondo-Heisenberg chain. In 
addition to the PDW phase discussed by Berg, Fradkin 
and Kivelson,^ we find that quantum frustration of the 
Heisenberg chain leads to the complex phase diagram 
shown in Fig. [TJ In the preceding sections we found that 







a. 

LL 






FIG. 1. Schematic phase diagram for the frustrated Kondo- 
Heisenberg chain. Here Ji and J2 respectively represent the 
nearest-neighbor and next-nearest neighbor exchange interac- 
tions and Jk the Kondo exchange interaction between the the 
spins of the frustrated Heisenberg chain and the ID electronic 
system. PDW denotes the pair-density-wave superconduct- 
ing phase, FSL the fractionalized spin Uquid (controlled by 
the Toulouse fixed point (TFP)), and D+LL the phase with a 
dimerized Heisenberg chain coexisting with a Luttinger liquid 
in the ID electronic system. The quantum phase transition 
between the D-I-LL phase and the PDW and FSL phases are 
most likely first order (see text). The phase boundary between 
the FSL and the PDW phases is of the Kosterlitz-Thouless 
class, for Jk weak enough, but may be first order for Jk large 
enough. We have not found any evidence of a phase transition 
at Jk — >■ 00. 

the phase diagram has the following phases: 

1. A PDW phase in the weak frustration regime, Ji ^ 
J2. This phase is controlled by a stable fixed point 
at large Jk and weak J2. 

2. For J2 > J2,c and Jk small there is a dimerized 
phase of the spin chain (with a spin gap) coexist- 



ing with a spin-1/2 Luttinger liquid in the ID elec- 
tronic system. The phase boundary between the 
dimerized spin chain-|-Luttinger liquid (D+LL) is 
most likely weakly first order since the system at 
Jk = and J2 = J2,c is at a fixed point with two 
marginally relevant operators. 

3. A fractionalized spin liquid (FSL) phase with a 
gapless (fractionalized) spin fluid coexisting with 
a decoupled charge sector of the ID electronic sys- 
tem. We investigated the nature of the quantum 
phase transition between the dimerized and the 
PDW phase and concluded that it most likely a 
first order transition. We also investigated the 
stability of the gapless spin fluid phase (charac- 
terized by a stable fixed point of the Toulouse 
form at finite Jk) and the nature of its correlation 
functions. We also showed that, for weak enough 
Jk, there is a Kosterlitz-Thouless (KT) type quan- 
tum phase transition from the gapless fractional- 
ized spin fluid to another phase which, due to a dan- 
gerous irrelevant operator, it is physically equiva- 
lent to the PDW phase of the weakly frustrated 
regime. This KT transition most likely becomes 
a first order transition at large enough Jk since 
in for Jk ^ 00 the system is unstable to a flow 
either towards the Toulouse fixed point (TFP) of 
towards the stable fixed point of the PDW phase. 
We have not found evidence a direct phase transi- 
tion between the dimerized phase and the gapless 
fractionalized spin fluid. This suggests that either 
the PDW phase "sneaks in" between this phase and 
the dimerized phase even in the strong frustration 
regime or, what is far more likely, that there is a 
direct first order phase transition. 

One of our motivations for looking at the role of frus- 
tration in the spin chain was to find out if a magnetic 
mechanism can give rise to an incommensurate PDW 
phase in a system without attractive interactions. As 
we noted in the introduction such a mechanism appears 
to be forbidden in ID systems with full SU(2) spin ro- 
tational invariance. Nevertheless such incommensurate 
phases are possible if the SU(2) spin symmetry is ex- 
plicitly broken down to its U(l) subgroup by a magnetic 
anisotropy interaction in the spin chain. 
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Appendix A: Solution at the Toulouse Point 



The Hamiltonian of Eqs. (|4.8p and (|4.9p can be treated 
using abelian bosonization as follows. We will first de- 
scribe the spin degrees of freedom of the electronic sys- 
tem in terms of a Bose field that we will denote by (p and 
whose dual field is -d. Their chiral components are, as 
usual, ipii^ = {ip + 'i9)/2 and ipL — {f ~ ^)/2 respectively. 
The chiral fields (pa^L are defined at the SU{2)i com- 
pactification radius, such that the chiral spin currents of 
the electronic system, Jg ^ and J^^s are given by the ex- 
pressions of the form of Eq. p.ll[) in terms of the chiral 
fields (fs (with s — R,L), 



^3,s 



1 

2'Ka 



j,=pi\/8iT-v3s 



1 



'3,s 



dx^s 



(Al) 



We will also need an expression for the chiral SU(2)2 
currents /" of Eq. (|4.6p in abelian bosonization. This is 
accomplished by combining two of the Majorana Fermi 
fields, S,i,s and ^2,3 into a single free Dirac (complex) 
Fermi field V's = {Q + i^l)/"^- The free Dirac fermion 
has a representation in terms of a chiral Bose fields (pn 
and (j)L 



1 



1 



fc = ;^(s+.^a ^,„ 



1 



1 



*' = 71<'^'^"^'' ■J2„ 



,iV47r0H 



-'iV4'7r0L 



(A2) 



In this language the chiral SU(2)2 currents of Ea. (l4.6p 
are given by 



7ra 
i 



^RFB.e 



±i\/4TTd 



--dx4>i 



If = ^^LFLe^'^^^^^ Il = ^dx^L (A3) 



where -Fi?,L are two (anticommuting) Klein operators. It 
is easy to see that the SU(2)2 chiral currents have scaling 
dimension 1 (as they should) and obey an SU(2)2 Kac- 
Moody current algebra. 

In this formulation the Hamiltonians "Hi and 'H2 of 



Eq. (l4.9p take the form 






H2 =Vi{dx^L? 



'-^ih^x^L 



V2{dx>PLT 



gill 

V2tt 



dx(f>RdxVL 



lnJlm 



^%Fr cos(V47r0fl 



8tt(pl) 
(A4) 



V2{dxVR)'^ 



9i\\ 
V27r 



dx<j>Ldx'PR 



27r-v/7ra 



^IFl cos(V47r(/)i 



Snipfi) 
(A5) 



The derivative couplings in the Hamiltonian Hi in 
Eq. (|A4p can be eliminated by means of a canonical trans- 
formation of the form 



(PR 



cosh a — sinh a 
— sinh a cosh a 



fii? 

$2L 



provided we choose 

tanh(2Q;) 



ffill 



7r-\/2(wi -I- V2) 



(A6) 



(A7) 



Similarly, the derivative couplings in the Hamiltonian "^2 
of Eq. (jA5|) can be eliminated by a canonical transforma- 
tion of the same form (with the same value of a) 



fR 



cosh a 
— sinh a 



— sinh a 
cosh a 



*1L 



(A8) 



For general values of gij| the cosine terms in Eqs. (jA4p 
and (|A5|) have non-trivial scaling dimension. The 
Toulouse point is defined as the value of gin for which 
the scaling dimension of the cosine operators is 1/2. This 
happens for 



tanh a^^ 



1 

71' 



9n 



An ^ 
y(«i 



V2) 



At the Toulouse point the Hamiltonian Hi 
Eq. (|4.9p becomes 

Hi + H2 =Ui [{dx^iRf + [dx^lLf] 
+U2 [{dx^2R)^ + (5,$2l)'] 

-'-Y&xeR-eLdxeL) 



(A9) 
H2 of 



„ ^- Cfl-Fi? cos(V47r$2L 
+ ^IFl cos(\/4^$2i?) 



(AlO) 



where ui = (2i;i — f2)/3 and U2 — (2u2 — vi)/i are the 
renormalized velocities at the Toulouse point. 

Since the cosine operators of Eq. (jA10| ) have scaling 
dimension 1/2, they can be refermionized (exactly as in 
the case of the two-channel Kondo problen>22i) . Let xi^) 
denote a Dirac Fermi field, whose chiral components are 



1 



XR 



V2 



(X2,_R + iXl,fl), XL 



1 



^/2 



(X2,L +iXl.L) 



(All) 
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respectively, where Xj,s are the two Majorana compo- 
nents (j — 1,2) with both chirahties (s = R,L) of the 
Dirac fermion. In turn, the chiral Dirac fermions XR and 
XL are related to the chiral bosons $2,i?, and ^2,l by an 
expression of the same form as in Ea. (jA2[) . i.e. 



This spectrum is massive with a single-particle gap of 
myl — (ui/w+) < m, where ui = (2wi — W2)/3 and 

V+ = {Vi +V2)/'S. 



1 



XR 



-.F, 



Re 



;-\/4¥<I>2 



1 



XL 



--Fj.e~ 



/47r$2 



(A12) 
where Fr and Fl are two Klein factors. 

In this basis, at the Toulouse point the Hamiltonian 
T-Ls of Eq. (|4.8p can be written as a sum of four decoupled 
Hamiltonians respectively consisting of the Luttinger liq- 
uid Hamiltonian for the Bose field $1 (and its dual field 
Gi), the Hamiltonians for the free Majorana fermions ^° 
and x^, and the Hamiltonian for the Majorana fermions 
^^ and X2 coupled to each other through a mass term, 

Us = HllI^i] + Hplf] + HfIx'] + HF[e,x'] (A13) 



where 



Ul 



nLLi^l] =y [{d.^lf 



(d.Qif 



HF[e 



«"1 /fO 



{eRd.eR - ad^a) 



(A14) 
(A15) 



tvi 



nF[x'] - ~ ^ {x'Rd.x'R - xi^.xi) (ai6) 



HFK^X^]=-^(c|5.C^-ei5.e 



^ [x^R^xX^R - xidxxl) 
im{eRxl-xUl) 



(A17) 



where we have defined the mass m = gij_/{2TTa). The 
Hamiltonian 'Hf[£,^,x'^] of Eq. (jA17P represents a theory 
of two free Majorana fermions coupled through the mass 
term shown in the last term on the right hand side of 
Eq. (IA17P . The spectrum of single-particle states of the 
Hamiltonian 'H_f[^'^,x^] is 



E±±ik) =±Mifc± 




(A18) 



Therefore, at the Toulouse point the Majorana 
fermions ^^ and x^ become massive and decouple from 
the other low energy degrees of freedom, represented 
by the massless boson i>i, and the massless Majorana 
fermions ^° and x^- 



The 2D classical Ising model and, equivalently, the ID 
(quantum) Ising model in a transverse field, is equiva- 
lent to a theory of Majorana fermions with a Majorana 
mass term which is tuned to zero at criticality.— i^ In 
the language of 2D conformal field theory, the Ising order 
parameter which we will denote by cr, and its Kramers- 
Wannier dual, the disorder operator which we will denote 
by fj,, respectively have non-vanishing expectation val- 
ues on each of the order and disorder phases of the Ising 
modelj^i^ In the fermionic language the two phases cor- 
respond to the two possible signs of the Majorana mass 
term (also known as the energy operator of the Ising 
model, which we will denote by e)^ In addition of label- 
ing the order and disorder phases, the order and disor- 
der operators twist the boundary conditions of the Ma- 
jorana fermion from periodic to anti-periodic, and hence 
are called twist-fieldSi-^- The scaling dimensions of the 
order and disorder operators at the Ising critical point, 
where the Majorana fermion is massless, is 1/8. In the 
construction outlined in this Appendix we encountered 
six Majorana fermions, ^", ^^, ^^, ^•^, x^ and x^- There- 
fore we will also have six associated twist fields a^, . . .a^ 
(in this order), and their associated Kramers- Wannier 
duals ^°,...,/i'^. These order and disorder twist fields 
enter in the definition of the observables of the frustrated 
Kondo-Heisenberg chain in this description. 
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